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Abstract 

Tt"1 ^ e s tudy a stochastic community model able to interpolate from a neutral regime to a niche partitioned regime upon varying a 

5s"^ single parameter, tuning the relative strength of intraspecific and interspecific competition. By means of a self-consistent approach, 

"; we obtain an analytical expression for the species abundance distribution, in excellent agreement with stochastic simulations 

O of the model. In the neutral limit, such calculation recovers the Fisher log-series. In the niche-dominated regime, the species 

^ abundance distribution develops a maximum at a finite population size, corresponding to the emergence of a carrying capacity. 

I We also investigate numerically species extinction-time distribution and compare it with the neutral case. We find that enhanced 

O competition increases the species lifetimes on average while suppressing their fluctuations. The results are discussed in view of the 

niche-neutral debate and of their potential relevance to field data. 
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1. Introduction 

Understanding the forces shaping ecological communi- 
ties and promoting species coexistence is a key longstanding 
issue in theoretical ecology. Niche theory (MacArthur and 
Levins, 1967a; Chesson, 2000; Chase and Leibold, 2003) 
and neutral theory (Caswell, 1976; Bell, 2001; Hubbell, 
2001) propose two alternative mechanisms for explaining 
long-term coexistence of species communities. 

Many observational studies have highlighted the impor- 
tance of specificities in resource exploitation and thus of 
niche differentiation for long-term coexistence (see Wright, 
2002; Silvertown, 2004, for two recent reviews). From a 
theoretical point of view, niche partitioning enhances com- 
petition with conspecific. As a consequence, species limit 
their own populations more than those of other species 
thus promoting coexistence (Chesson, 2000). A simple 
mathematical formalization of this idea is provided by the 
Lotka-Volterra competition model (MacArthur and Levins, 
1967a), where coexistence emerges whenever intraspecific 
interactions overtake interspecific ones. 

Apparently on the other side of the Moon, ecological neu- 
tral theory (Hubbell, 2001) completely disregards species 
differences and assumes functional equivalence at the in- 
dividual level. Within this theoretical approach, stochas- 
ticity is the leading ecological force, and species coexis- 
tence is the result of a dynamical balance between immi- 



gration/speciation processes and extinction. Thanks to the 
drastic assumptions of neutral theory, ecological patterns 
such as species abundance relations can be analytically 
computed in this framework. These predictions, depending 
on few parameters, fit surprisingly well observational data 
from tropical forests (see, e.g., Hubbell, 2001; Volkov et al., 
2003, 2005). 

Due to its strong departure from traditional theoretical 
approaches, the neutral theory elicited a heated debate in 
community ecology about its validity and interpretation. A 
complication is that, in many cases, niche-based and neu- 
tral models yield compatible fits of biodiversity patterns 
(see, e.g., Chave et al., 2002; McGill, 2003; Mouquet and 
Loreau, 2003; Tilman, 2004; McGill et al., 2006). 

Here, we do not enter this debate, but embrace the view 
that "niche and neutral models are in reality two ends of a 
continuum with the truth most likely in the middle" (Chase 
and Leibold, 2003). Indeed, the ecological forces underly- 
ing niche and neutral models are not by no means mutu- 
ally exclusive. Following a similar view, different authors 
have been proposing models synthesizing niche-based and 
neutral mechanisms (Tilman, 2004; Leibold and McPeek, 
2006; Gravel et al., 2006; Adler et al, 2007; Kadmon, 2007; 
Pigolotti and Cencini, 2010; Haegeman and Loreau, 2011; 
Noble et al., 201 la, b). A desired property of these models 
is the possibility to recover pure neutrality as a limiting 
case (Alonso et al., 2006), allowing for studying the transi- 
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tion from niche-dominated to neutral regime when one (or 
more) ecological parameter is varied. A related issue con- 
cerns the robustness of the neutral theory. In particular, a 
general understanding of whether small non- neutral effects 
can undermine neutral predictions is still lacking. Finally, 
niche-neutral models could help our understanding of the 
dynamical features and long time behavior of ecosystems. 
Indeed, while neutral theory in principle allows for extrap- 
olating ecosystem behavior at large time-scales (Pigolotti 
et al, 2005; Azaele et al, 2006; Bertuzzo et al., 2011), an 
open issue is that predicted species lifetimes tend to be too 
short compared to fossil records estimates (see, e.g. Lande 
et al, 2003; Ricklefs, 2003; Nee, 2005; Ricklefs, 2006; Allen 
and Savage, 2007) 

A difficulty with the above program is that usually such 
mixed models are very hard to analyze mathematically (see, 
e.g., Tilman, 2004). In this respect, surely, one of the rea- 
sons of success of the neutral theory has been its mathe- 
matical tractability, i.e. the fact that closed expression for 
distribution such as species abundance distributions (SAD) 
can be obtained and easily compared with data (Volkov 
et al., 2003, 2005). This is usually not possible in models 
where the neutral hypothesis is broken. Moreover, break- 
ing the neutral hypothesis often leads to a proliferation of 
parameters, making general comparisons problematic. 

As discussed by Adler et al. (2007), neutrality can be vi- 
olated both at the level of fitness inequalities among the 
species and/or by the presence of stabilizing processes caus- 
ing species to limit themselves more than they limit others. 
These are the two main mechanisms for coexistence identi- 
fied by Chesson (2000). As for the former kind of violation, 
Zhang and Lin (1997); Zhou and Zhang (2008) have shown 
that even small differences in fitness in an otherwise neutral 
model cause strong reduction of biodiversity. This result is 
not unexpected as fitness differences, in the absence of sta- 
bilizing effects, clearly lead to competitive exclusion. As for 
the second kind of violation of neutrality, quite interesting is 
the approach followed by Haegeman and Loreau (2011) who 
added demographic stochasticity and immigration to clas- 
sical competitive Lotka-Volterra dynamics with symmetric 
interactions. A slightly more general model, developed in 
(Noble et al., 2011a, b; Noble and Fagan, 2011c), incorpo- 
rates both fitness inequalities and stabilizing mechanisms. 
In a nutshell, Noble et al. added a frequency-dependent 
birth rate to the Moran model in a such a way that it re- 
duces to Hubbell's neutral model in an appropriate limit 
and recovers standard Lotka-Volterra phenomenology in 
the deterministic limit. They mostly explored ecosystems 
with few species, aiming at understanding niche stabiliza- 
tion effects at intermediate time scales. 

In this paper we investigate a simple variant of the model 
introduced by Noble et al. (2011b) (see also Noble and Fa- 
gan, 2011c). In particular, we focus on communities sub- 
ject to speciation, death and reproduction where fitness 
inequalities are disregarded, while stabilizing mechanisms 
are retained via symmetric interactions as in Haegeman 
and Loreau (2011). Indeed field data studies (see. e.g., 



Volkov et al., 2009, for a recent study in this direction) 
support the idea that the most relevant deviation from 
neutrality consists in differences between intraspecific and 
interspecific competition. Such differences are usually ex- 
plained in terms of species-specific resource consumption 
or more sophisticated mechanisms like Janzen-Connell ef- 
fects (Janzen, 1970; Connell, 1971). 

By means of a self-consistent computation, we derive a 
closed expression for SAD, depending on a single parameter 
which quantifies the deviation from neutrality. In the neu- 
tral limit, the classical Fisher log-series Hubbcll (2001) is 
recovered. Out of the neutral limit, competitive interactions 
are essentially incorporated into a single self-consistently 
derived quantity, which modifies the Hubbell's fundamen- 
tal biodiversity parameter. Predicted SAD are successfully 
compared with numerical simulations. Moreover, we study 
the species extinction time statistics as a function of the 
strength of the stabilizing effects. 

Model details are summarized in Sec. 2. Section 3 presents 
the main results: the analytical derivation of the SAD with 
numerical validation, and a numerical study of the species 
extinction times distribution as a function of the deviation 
from neutrality. Section 4 is devoted to discussing the re- 
sults. The Appendix contains some technical details of the 
self-consistent derivation presented in Sec. 3. 

2. Model 

For the sake of simplicity, we start by introducing the 
neutral version of the model, which is the standard Moran 
model for a community with a fixed number, J, of indi- 
viduals. At each timestep an individual, randomly chosen 
among the J composing the community, is killed. With 
probability v, it is replaced with an individual belonging 
to a new species, not present in the system (speciation). 
With probability (1 — v), it is replaced with a copy of one 
of the J — 1 individuals already present in the community 
(reproduction). Denoting with rii the number of individu- 
als of species i — 1,...,S (S being the number of species 
currently in the system) present at time t, the death of an 
individual of species j followed by the birth of one of species 
i happens at rate 



(!-")■ 



J J -I ' 



(1) 



with ni,n,j > 1 . 

We now move to the general non-neutral case. Following 
Noble et al. (2011b), non- neutrality is introduced by bias- 
ing reproduction. We introduce the frequency-dependent 
weights 

w t = r 4 e~^o=i Ql3nj , (2) 

and substitute the term m/ (J — I) of Eq. (1) with 

u>i rii ,. 

^s ■ ( 3 ) 

where r^ models the intrinsic fitness advantage of species i 
and ciij the competition among the species present in the 



community. With equal fitness r, — r and uniform competi- 
tion ciij = aj , Eq. (3) reduces to the neutral value nj/( J— 1) 
(by definition X^fe=i Hfc = J— 1 as the dead individual docs 
not participate to reproduction). In a different limit, one 
can draw a correspondence with the Lotka-Volterra com- 
petitive model by expanding (3) for small a^ and neglecting 
fluctuations (Noble et al., 2011b; Noble and Fagan, 2011c). 
Specifically, we consider a simplified variant of the model 
by imposing equal fitness r, = r while retaining niche- 
stabilization (Adler et al., 2007) through a fully symmetric 
competition matrix a%j = aSij + b(l — 5ij), as in Haegcman 
and Loreau (2011). With these simplifications, Eq. (3) re- 
duces to 
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where we used X^i n j = (J ~ 1) — n i- Parameters a and 
b tune the weights of intraspecific and interspecific com- 
petition, respectively. Finally, by defining c = a — b, the 
frequency-dependent rates can be written as 
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where c fully specifies the deviation from the neutral case, 
corresponding to c = 0. 

We focus on the case c > 0, meaning intraspecific com- 
petition being stronger than interspecific one (i.e. a > b), 
as observed in real data (Volkov et al., 2009). 

3. Results 

In this section we present the main results. First, we 
show that the species abundance distribution (SAD) of the 
model can be analytically computed by means of an ap- 
proach similar to that used in Volkov et al. (2003) for the 
neutral model, complemented with a self-consistent ansatz 
for the treatment of non-neutral competitive interactions. 
Second, we study the dynamical properties of the model 
by numerically investigating the single species extinction 
times statistics. 

3.1. Species abundance distribution 

Within the model introduced in Sect. 2, the dynamics of 
the population of a species, labeled i, is ruled by the master 
equation 

d 



dt 
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The latter, instead, explicitly accounts for competition with 
the other species and reads 
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The complication with respect to the neutral case is that, 
now, the master equation for the focal species depends on 
the populations of all the others. Notice that, in Eq. (8), the 
birth probability a/ 4 '(n) is decreased by the factor (1 — v) 
to account for speciation, which from the point of view of 
the focal species corresponds to the possibility to give birth 
to an individual of a new species. Following Volkov et al. 
(2003) we set &q = v (accounting for the appearing of the 
species) and b_\ = 0. It is worth underlining that, with 
this representation, the total number of individuals is not 
fixed, which is a good approximation when both the total 
population size J and the number of species S are large. 

Imposing the detailed balance, the equilibrium distribu- 
tion can be written as 
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However, unlike the neutral case, expression (9) does not 
solve the problem as it represents just a formal solution, 
valid for given values of the populations of the other species 
m (with I =/= i), which are unknown. The normalization 
factor 8 is introduced in analogy with the "fundamental 
biodiversity number" in Hubbell's model (Hubbcll, 2001) 
and will be made explicit below. 

To solve the problem, we can now exploit the species- 
level symmetry of the model and make the ansatz that all 
species have the same equilibrium distribution, thus drop- 

(i) 

ping the dependence on i, i.e. p k = pf.. Further, if parame- 
ters are such that S is sufficiently large, we can assume that 
the ni are uncorrclatcd random variables. These two as- 
sumptions are sufficient to find the solution to our problem 
self-consistcntly. Upon defining the generating function of 
the equilibrium distribution 
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and under the two assumptions above, we can simplify the 
denominator of (9), yielding 



Pk = 



J k (l-v) k e-^ k - 1 ^ 



k[-Sp'{c)] k 
where the prime denotes derivative, so that p'(c 
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where Pn (t) is the probability for species % of having n % = J2k=iPkke c ■ Notice that Sp (0) — S(n) — J, so that 



n individuals at time t. The death and birth probabilities 
d„ and bjl depend on the species populations and can be 
expressed as follows. The former is completely neutral and 
determined by the relative population of species i 

4 J) = n -j ■ (7) 



we can rewrite the equilibrium distribution more com 
pactly as 
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Pk = 
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where g = p'(0)/p'(c). Equation (12) provides an analytical 
expression for the {pk}'s as a function of the deviation 




Fig. 1. (A) Species abundance distribution at varying the competition intensity c as in the legend. Inset: zoom of the same curves in log-lin 
scale to highlight the emergence of a mode. Symbols refer to simulations and solid curves to the prediction (12). (B) The (c,v) parameter 
plane. The grey region corresponds to parameter values for which the condition (15) is fulfilled, i.e. where a peak in the SAD is present. The 
symbols correspond to parameters as in the simulations shown in (A) but for the neutral value c = here not represented. 



from neutrality c, and constitutes our central result. In 
order to interpret p k as the SAD we have to impose the 
normalization condition ^2 k=1 kpk — J ■ As shown in the 
Appendix, this normalization leads to the following exact 
expression for the biodiversity number 

e = ^- v d3) 

s(i-^) 

In the neutral limit (c = 0), we have g — 1 so that 8 is 
equal to the neutral value Jvj(\ — v) and the SAD (12) 
reduces to the Fisher log-series (Hubbell, 2001; Volkov et 
al., 2003). Conversely, in the general case c ^ 0, g > 1 
and the biodiversity number is enhanced by with respect 
to the neutral case due to competition. In the Appendix, 
we show that the unknown quantity g can be obtained as 
the solution of the transcendental equation 
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In Figure 1A we show SADs from the model at v = 10 -3 
and different values of c as indicated by symbols in the c, v 
plane shown in Fig. IB, with the addition of the neutral 
value c = not represented therein. In all cases, numerical 
data (symbols) are in excellent agreement with the analyti- 
cal prediction (12) represented in solid lines, validating the 
self-consistent approach. 

For small values of c, the SAD is monotonically decreas- 
ing with n, similarly to the (neutral) Fisher log-series but 
for faster (Gaussian) fall-off for large populations. Interest- 
ingly, for larger values of c, SAD curves develop a maximum 
at a finite number of individuals, as emphasized in the in- 
set of Fig. 1A. A study of the function (12) reveals that the 
condition for developing a maximum at n > 1 is 
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corresponding to the shaded area in Fig. IB. 

We now discuss how the average number of species (S) = 
J2kP k depends on the deviation from neutrality. As this 



number is proportional to the number of individuals J, in 
Fig. 2A we plot the fraction (S) / J as a function of the spe- 
ciation rate v and the intensity of competition c. For any 
fixed value of the speciation rate z/, the average number of 
species grows with c, meaning that the system is able to 
sustain a larger diversity with respect to the neutral case 
thanks to the stabilizing effect of niches. In Fig. 2B, we 
compare the prediction obtained by summing the terms 
of Eq. (12) with numerical simulation of the model, again 
showing an excellent agreement between theoretical pre- 
diction and simulation of the stochastic model. 

3.2. Species extinction times distribution 

We now focus on dynamical, rather than static, patterns 
which provide useful information on the differences between 
neutral and non- neutral theories (Pigolotti et al, 2005; Aza- 
ele et al, 2006; Allen and Savage, 2007; Bertuzzo et al., 
2011). In particular, it is interesting to compare features of 
the dynamical balance between species introduction, and 
extinction as a function of the competition strength c. As 
in the model species originate from a single individual, we 
are interested in the statistics of the time it takes for a 
new species (introduced via a single individual) to become 
extinct. The results are presented in unit of generations, 
meaning that t = 1 corresponds to J iteration steps of the 
model of Sec. 2. 

In the neutral case c = 0, by applying ideas from the 
theory of branching processes (Harris , 1989), it is possible 
to derive an analytical expression for the probability of a 
species to survive a time t (from the time of its introduction 
in the system), which is given by (Pigolotti et al, 2005) 

For v <C 1) such distribution displays a power- law tail 
£~ 2 , followed by an exponential cutoff set by v. Notice 
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Fig. 2. (A) Average number of species normalized by the community size, (S)/J, in the (c, v) parameter space, as predicted by summing the 
terms of Eq. (12). (B) Same quantity as a function of v for different values of c, increasing in the arrow direction. Symbols are results from 
simulations, in excellent agreement with the prediction in red lines. Grey lines correspond to values of c different from those of the simulations. 



that the power-law t~ 2 can be modified by dispersal 
properties different from the global one here considered 
(Bertuzzo et al., 2011). For c > 0, no closed expression for 
the extinction-time distribution is available. Indeed, also 
a simple quadratic nonlinearity in the birth (or death) 
rates for a single species makes the problem analytically 
intractable, but for clever approximations as discussed in 
a series of works (see, e.g. Norden, 2001; Newman et al., 
2004; Doering et al, 2008; Parsons et al, 2008). For this 
reason, we limit ourself to a numerical study. In Figure 3, 
we show the extinction-time distribution, p e (£), for differ- 
ent values of c fixing v. Increasing c has a non-trivial effect 
on the shape of the extinction-time distribution. 

Qualitatively, one can identify three regimes. At short 
times, the distribution is only weakly dependent on c, i.e. 
by the presence of competition, and essentially reproduces 
the neutral result. At intermediate times, the extinction- 
time probability increases with c. Finally, at large extinc- 
tion times it decreases with c by developing an exponential 
cutoff much steeper than predicted by formula (16). Fitting 
the exponential rate of decrease suggests a dependency of 
extinction probability of the formp e (i) ~ c~ Qt , at least for 
small c (inset of Fig. 3). 

As suggested by niche-stabilization arguments, the net 

effect of these different regimes on the average extinction 

time is always positive: the average lifetime of a species in 

our model is a growing function of c (see Fig. 4) , which can 

exceed up to 2 or 3 order of magnitude (depending on v) 

the neutral value 

ln(l/i/) 
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derived from Eq. (16). It is interesting to notice that, while 
enhancing the intensity competition (larger c) increases the 
average extinction time, it decreases its relative fluctua- 
tions as confirmed by the behavior of the coefficient of vari- 
ation, CV = y/(t 2 ) - (i) 2 /(t), shown in the inset of Fig. 4. 
This means that, while in the neutral case when v <C 1 ex- 
tinction times can vary greatly around the average value, 




Fig. 3. Probability density function, p e (t), of the extinction times 
for J = 10 4 , v = 5 ■ 10~ 4 and different values of the competition 
parameter c, as in the legend. For c = the prediction of (16) is 
perfectly verified. Inset: behavior of the exponent of the exponential 
tail of pdf as a function of c. The solid line displays the fit a ln(e) + b 
with a ss 0.000312 and b = 0.000455. 

niche stabilization has also the effect of making the aver- 
age time a better prediction for the lifetime of a randomly 
chosen species. 



4. Discussions 

We studied a community model incorporating neutral 
demographic stochasticity with niche stabilization. The 
model belongs to a general class of niche-neutral mod- 
els discussed in the recent literature (Gravel et al., 2006; 
Adler et al., 2007) and, in particular, is built upon two spe- 
cific models originally proposed by Haegeman and Loreau 
(2011); Noble et al. (2011b); Noble and Fagan (2011c). 
We assumed that violation of neutrality is controlled by 




Fig. 4. Average extinction time as a function v for different values 
of c, as in the legend. For c = 0, the prediction of Eq. 17 perfectly 
agrees with the data. Inset: coefficient of variation of the extinction 
time, defined as CV = ({t 2 } — (i) 2 ) 1 / 2 /(£), as a function of v for the 
same values of c. 

a single parameter, c, weighting the relative intensity of 
intra- versus interspecific competition. We focused on the 
case in which effects of interspecific interactions are weak 
compared with intraspecific ones (i.e. c > in our model), 
a feature well supported by data (see, e.g., Condit et 
al, 1992; Harms ct al, 2000; Volkov et al., 2009) and by 
mechanisms such as the well known Janzen-Conncll effects 
(Janzcn, 1970; Conncll, 1971). For the sake of simplicity, 
as in (Haegeman and Loreau, 2011), we assumed that, 
even in non-neutral cases, species are symmetric and with 
no intrinsic fitness advantages. The model depends also 
on rate v of introduction of new species by speciation or 
immigration. 

We proposed a self-consistent calculation which al- 
lowed us to straightforwardly derive an analytical ex- 
pression for the species abundance distribution, Eq. (12). 
Self-consistency avoids some of the problems of other ap- 
proaches like pcrturbative techniques, where one has to 
figure out the relative magnitude of the two small parame- 
ters v and c. While we focused here on one specific source 
of non-neutrality, one can expect that other non-neutral 
(symmetric) models can be solved using similar ideas. 
Indeed, species permutation symmetry is the crucial as- 
sumption to compute the SADs, rather than the specific 
way niches are introduced. However, we remark that even 
in the presence of species permutation symmetry, the self- 
consistent ansatz may not lead to the correct result when 
this symmetry is spontaneously broken (Borile et al, 2012). 

Interestingly, in a large portion of the parameter space 
(c, v), the relative abundance distribution curves show a 
maximum for intermediate abundance classes, similarly 
to many instances of SAD curves measured on the field 
(Hubbcll, 2001). Specifically, such mode appears when in- 
traspecific competition parameter c exceeds a critical value 
(predicted by Eq. (15)) depending on the rate v of intro- 



duction of new species. We remind that spatially implicit 
implementations of the neutral model, which include some 
level of dispersal limitation, have been shown to provide 
very good fit of such field-measured SAD curves (Hubbell, 
2001; Volkov et al., 2003). Without dispersal limitation, 
good fits have been obtained by combining neutral com- 
munity models with specific forms of density-dependent 
reproduction rates (Volkov et al., 2005). 

Within the model here investigated, the presence of a 
peak in SAD curves corresponds to the natural emergence 
of an effective carrying capacity for each niche, so that the 
population of each species fluctuates around a well defined 
average value. It should be remarked, however, that the tail 
of the SAD distribution, i.e. the probability of encountering 
a very abundant species, falls off in a much sharper way for 
large n than in Hub-bell's spatially-implicit neutral model, 
due to the presence of a term proportional to exp(— en 2 ) in 
Eq. (12). Typical datasets well fitted by the neutral theory 
(Hubbcll, 2001; Volkov et al., 2003) support the presence of 
a fatter log-normal-like tail. As a consequence, best fits of 
those datasets with formula (12) (not shown) are biased by 
such long tails and are thus realized close or at the neutral 
limit (i.e. for c — > 0), where the peak disappears. 

The emergence of this Gaussian cutoff, conflicting with 
classic datasets, is likely due to the fact that the niche- 
stabilization mechanisms incorporated in our model are 
particularly strong in the absence of dispersal limitation. It 
is indeed reasonable to expect that the sharp fall-off should 
be a quite general feature for global dispersal models in 
which an effective carrying capacity for any single species 
emerges and should not depend too much on model de- 
tails. We conjecture that including some form of dispersal 
limitation with a competition acting on a finite range may 
lead to fatter tails. Similar features were observed in Chave 
et al. (2002) while comparing global and limited dispersal 
models with density dependence or other stabilizing mech- 
anisms, such as tradeoffs. In this perspective, it would be 
interesting to extend the model to incorporate, even in a 
spatially implicit form, some degree of dispersal limitation. 
Unfortunately, this extension makes the analytical treat- 
ment much harder and it is left for future investigations. 
Another way fatter tails may arise would be to break the 
symmetric hypothesis and letting different species having 
different carrying capacity. 

Studying dynamic patterns of species lifetimes (Sect. 3.2) 
revealed that the main effect of violating neutrality is to 
suppress the large fluctuations of lifetimes, typical of neu- 
tral dynamics (Pigolotti et al, 2005). On the one hand, en- 
hanced self-competition prevents the possibility of a species 
to achieve a very large population size and thus lifetimes 
much longer than the average, which are exponentially sup- 
pressed. Similar exponential distributions of lifetimes have 
been observed in fossil data, a fact often explained in terms 
of the Red Queen effect in a changing environment (Van 
Valen, 1973). On the other hand, a species is favored when 
rare making it less prone to extinction by demographic 
stochasticity when its population is small. These effects arc 



"weighted by the competition parameter, as a result the av- 
erage species lifetime results increased with respect to the 
neutral expectation (17). 

As observed by many authors (see, e.g. Lande et al, 2003; 
Ricklefs, 2003; Nee, 2005; Ricklefs, 2006; Allen and Savage, 
2007), one of the problem with the neutral theory relies 
on the fact that average species lifetimes are typically too 
short, up to some order of magnitude, compared to those 
estimated from fossil records. The origin of such discrep- 
ancy is likely due to the point speciation mode, typically 
implemented in neutral models (Hubbell, 2003). As sug- 
gested by Allen and Savage (2007), a possible way-out for 
overcoming this problem is to allow for larger incipient- 
species abundances. Moreover, the same authors have also 
shown that allowing for some environmental stochasticity 
can decrease the lifetime of very abundant species, which 
is another issue with the neutral prediction for the species 
lifetimes distribution. 

Furthermore, our model shows that niche-stabilization 
induced by intraspecific interactions larger than interspe- 
cific ones, for a given value of v, both increases the average 
lifetime (up to 2 or 3 order of magnitudes depending on v 
and c) and suppresses large fluctuations. In particular, the 
smaller the value of v, the larger is the effect of competition. 
Thus, competition allows the system to sustain a larger di- 
versity (Fig. 2) of longer-lived (Fig. 4) species with respect 
to a purely neutral community and, at the same time, sup- 
pressing the probability of species with exceedingly large 
lifespan with respect to the average. 
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Appendix A. Computing g and 6 

Here, we derive the expression for the normalization 9 
presented in Eq. (13), and discuss a semi-analytical method 
to estimate the parameter g, which is self-consistently de- 
fined as 

■ p'(c) Ek=i( l -v) k g k e- §Hk+1) 

For the sake of notation simplicity we denote with Af and 
T>, respectively, the numerator and denominator in the right 
hand side of expression (A.l). By rearranging the indices 
in the summations, one can show that they satisfy 



Z> = £(l-i/) 



k k -ck(k+l)/2 _ 

y e 



1 -v 



(A.3) 



V 



N-g{l-v) 



.9(1 - v) 
Substituting (A. 2) into (A.l) yields an expression for V 



(A.2) 



fe=i 



The normalization 6 can be then derived by imposing 
the condition J^ fc kp(k) = J. By direct substitution, one 
obtains J2k kp(k) = ON. Then, using the fact that N = 
gV = g(l — v)jv yields the result (13). 

Let us now discuss how estimating g in the relevant pa- 
rameter range of c ^C 1 and v <C 1. In such limit, one can 
approximate very closely the series in (A.3) with an integral 



1 -v 



crfc 



/oo r~ 

dk (1 - v) k g k e-^ k{k+l) = J— (A.4) 

ex P{ l^Mi^M! } 



3c - 2 ln(.g(l - u)) \ _ _ f [c - 2 ln((l - v)g)] 2 



2y/2c 



8c 



Equating the first and the last term in the above expression 
leads to Eq. (14), that can be solved for g semi- analytically 
by standard methods. The dependence of g on v and c is 
shown in Fig. A.l. 




Fig. A.l. Contour plot of the quantity g in the (c, y)-plane, as com- 
puted by solving Eq. (14) semi-analytically. 

We conclude by discussing possible analytical expres- 
sions for g. In the neutral limit of c = 0, one has g = 1, it is 
thus tempting to build up an expansion to obtain a closed 
expression for g at least in the limit c«l. However, such 
expansion in Eq. (A.4) requires additional assumptions on 
the relative magnitude of the two small parameters c and 
v. For example, the argument of the error function is very 
small when v -C c -C 1, but very large when c <C v <C 1, 
so that in the first case a Taylor expansion is appropri- 
ate, while in the second one has to perform an asymptotic 
expansion. More formally, one can show that the function 
g{v,c) is not analytic at the point v = c = 0, so that one 
can not perform a Taylor expansion in the two small pa- 
rameters. This also suggests that, more in general, other 
near-neutral models can be hard (or impossible) to treat 
with perturbative methods due to the interplay between 
the two small parameters. 



References 

Adler, P., J. Hillcrislambers, and J. Levine. 2007. A niche 
for neutrality. Ecol. Lett. 10, 95-104. 

Alonso, D., R. S. Etienne, and A. J. McKane. 2006. The 
merits of neutral theory. Trends Ecol. Evol. 21, 451-457. 

Allen, A. P. and V.M. Savage. 2007. Setting the absolute 
tempo of biodiversity dynamics. Ecol. Lett. 10, 637646. 

Azaele, S., S. Pigolotti, J. R. Banavar, and A. Maritan. 
2006. Dynamical evolution of ecosystems. Nature 444, 
926-928. 

Bell, G. 2001. Neutral Macroecology. Science 293, 2413- 
2418 

Bcrtuzzo, E., S. Suweis, L. Mari, A. Maritan, I. Rodriguez- 
Iturbe and A. Rinaldo. 2011. Spatial effects on species 
persistence and implications for biodiversity. Proc. Natl. 
Acad. Sci. 108, 4346-4351. 

Borilc, C, M. A. Munoz, S. Azaele, J. Banavar, and A. Mar- 
itan. 2012. Spontaneously Broken Neutral Symmetry in 
an Ecological System. Phys. Rev. Lett. 109, 038102. 

Caswell, H. 1976. Community structure: a neutral model 
analysis. Ecol. Monogr. 46, 327-354. 

Chase, J. M., and M. A. Leibold. 2003. Ecological Niches: 
Linking Classical and Contemporary Approaches. The 
University of Chicago Press. 

Chave, J., H. C. Muller-Landau, and S. A. Levin. 2002. 
Comparing Classical Community Models: Theoretical 
Consequences for Patterns of Diversity. Am. Nat. 159, 
1-23. 

Chcsson P. 2000. Mechanisms of manteinance of species 
diversity. Annu. Rev. Ecol. System. 31, 343-66. 

Condit, R., S. P. Hubbell, and R. B. Foster. 1992. Recruit- 
ment near conspecific adults and the maintenance of tree 
and shrub diversity in a Neotropical forest. Am. Nat. 140, 
261-286. 

Conncll, J. H. 1971. On the roles of natural enemies in pre- 
venting competitive exclusion in some marine animals 
and in rain forest trees. Pages 298-312 in P. J. den Boer 
and G. R. Gradwell, eds. Dynamics of populations. Pro- 
ceedings of the Advanced Study Institute on "Dynamics 
of Numbers in Populations," Oosterbeek, Netherlands. 
Pudoc, Wageningcn. 

Docring C. R., K. V. Sargsyan, and L. M. Sanders. 2008. 
Extinction times for birth-death processes: exact results, 
continuum asymptotics, and the failure of the Fokkcr- 
Planck approximation. Mutliscale Model. Simul. 3, 283- 
299. 

Gravel, D., C. D. Canham, M. Beaudet, C. Messier. 2006. 
Reconciling niche and neutrality: the continuum hypoth- 
esis. Ecol. Lett. 9, 399-409. 

Haegeman, B., and M. Loreau. 2011. A mathematical syn- 
thesis of niche and neutral theories in community ecol- 
ogy. J. Theor. Biol. 269, 150-165. 

Harms, K. E., S. J. Wright, O. Caldcron, A. Hernandez, 
and E. A. Hcrre. 2000. Pervasive density-dependent re- 
cruitment enhances seedling diversity in a tropical forest. 



Nature 404, 493-495. 

Harris, T. E. 1989. The Theory of Branching Processes. 
Dover Univ. Press, New York. 

Hubbell, S. P. 2001. The Unified Neutral Theory of Bio- 
diversity and Biogeografy. Princeton University Press, 
Princeton, NJ. 

Hubbell, S. P. 2003. Modes of speciation and the lifespans 
of species under neutrality: a response to the comment 
of Robert E. Ricklefs. Oikos 100, 193-199. 

Janzen, D. H. 1970. Herbivores and the number of tree 
species in tropical forests. Am. Nat. 104, 501-528. 

Kadmon, R., and Allouche O. 2007. Integrating the effects 
of area, isolation, and habitat heterogeneity on species 
diversity: a unification of island biogeography and niche 
theory. Am. Nat. 170, 443-454. 

Landc, R., S. Engen and B. E. Saether, 2003. Stochastic 
Population Dynamics in Ecology and Conservation. Ox- 
ford University Press, Oxford, UK. 

Leibold, M. A., M. A. McPcck. 2006. Coexistence of the 
niche and neutral perspectives in community ecology. 
Ecology 87, 1399-1410. 

MacArthur, R., and R. Levins. 1967. The limiting simi- 
larity, convergence, and divergence of coexisting species. 
Am. Nat. 101, 377-385. 

McGill, B. J. 2003. Strong and weak tests of macroecolog- 
ical theory. Oikos. 102:679-685. 

McGill, BJ., BA. Maurcr, M. D. Weiser. 2006. Empirical 
evaluation of neutral theory. Ecology 87, 1411-1423. 

Mouquet, N., and M. Loreau. 2003. Community Patterns in 
SourccSink Mctacommunitics. Am. Nat. 162, 544-557. 

Nee, S. 2005. The neutral theory of biodiversity: do the 
numbers add up? Funct. Ecol. 19, 173-176. 

Newman T. J., J.-B. Ferdy, and C. Quince. 2004. Extinc- 
tion times and moment closure in the stochastic logistic 
process. Theor. Pop. Bio. 65, 115-126. 

Noble, A. E., N. M. Tcmme, W. F. Fagan, and T. H. Keitt. 
2011. A sampling theory for asymmetric communities. J. 
Theor. Biol. 273, 1-14. 

Noble, A. E., A. Hastings, and W. F. Fagan. 2011. Multi- 
variate Moran process with Lotka-Volterra phenomenol- 
ogy. Phys. Rev. Lett. 107, 228101. 

Noble, A. E., and W. F. Fagan. 2011. A unification of niche 
and neutral theories quantifies the impact of competition 
on extinction. arXiv preprint arXiv:1102.0052. 

Norden R. H. 2001. On the distribution of the time to 
extinction in stochastic logistic population model. Adv. 
Appl. Probab. 14, 687-708. 

Parsons T. L., C. Quince, and J. B. Plotkin. 2008. Ab- 
sorption and fixation times for neutral and quasi-neutral 
populations with density dependence. Theor. Pop. Biol. 
74, 302-310. 

Pigolotti, S., A. Flammini, M. Marsili, and A. Maritan. 
2005. Species lifetimes for simple models of ecologies. 
Proceedings of the National Academy of Science 102, 
15747-15751. 

Pigolotti, S., and M. Cencini. 2010. Coexistence and inva- 
sibility in a two-species competition model with habitat 



preference. J. Theor. Biol. 265, 609-617. 

Ricklefs, R. E. 2003. A comment on Hubbell's zero-sum 
ecological drift model. Oikos 100, 185-192. 

Ricklefs, R.E. 2006. The unified neutral theory of biodiver- 
sity: do the numbers add up? Ecology, 87, 1424-1431. 

Silvcrtown J. 2004. Plant coexistence and the niche. Trends 
Ecol. Evol. 19, 605-611. 

Tilman, D. 2004. Niche tradeoffs, neutrality, and commu- 
nity structure: A stochastic theory of resource compe- 
tition, invasion, and community assembly Proc. Natl. 
Acad. Sci. 101, 10854-10861. 

Van Valen, L. 1973. A new evolutionary law. Evol. Theory 
1, 1-30. 

Volkov, I., J. R. Banavar, S. P. Hubbcll, and A. Maritan. 
2003. Neutral theory and relative species abundance in 
ecology. Nature 424, 1035-1037. 

Volkov, I., J. R. Banavar, F. He, S. P. Hubbell, and A. Mar- 
itan Volkov, Igor, Jayanth R. Banavar. 2005. Density de- 
pendence explains tree species abundance and diversity 
in tropical forests. Nature 438, 658-661. 

Volkov, I., J. R. Banavar, S. P. Hubbell, and A. Maritan. 
2009. Inferring species interactions in tropical forests. 
Proc. Natl. Acad. Sci. 106, 13854-13859. 

Wright, S. J. 2002. Plant diversity in tropical forests: a 
review of mechanisms of species coexistence. Oecologia 
130, 1-14. 

Zhang, D.-Y., and K. Lin. 1997. The effects of competi- 
tive asymmetry on the rate of competitive displacement: 
how robust is Hubbells community drift model? J. Theor. 
Biol. 188, 361-367. 

Zhou, S. R., and D. Y. Zhang. 2008. A nearly neutral model 
of biodiversity. Ecology 89, 248-258. 



